      IO       1.0         MARCH     1995      OPEN FILES VERS, MESH, INCON, GENER, SAVE, LINEQ, TABLE, XAMD, XAMDF, PLOT, TRACE

      TOUGH-LGM1.1      31 OCTOBER   1993      MAIN PROGRAM FROM TOUGH2 VERSION 1.0 (29 MARCH 1991)

      INPUT    1.0         MARCH     1995      READ ALL DATA PROVIDED THROUGH FILE *INPUT*
      FLOP     1.0      11 APRIL     1991      CALCULATE NUMBER OF SIGNIFICANT DIGITS FOR FLOATING POINT ARITHMETIC
      RFILE    1.2         MARCH     1995      INITIALIZE DATA FROM FILES *MESH* OR *MINC*, *GENER*, *INCON* AND *XAMD*
      INDATA   1.2      29 SEPTEMBER 1993      PROVIDE PRINTOUT OF MOST DATA PROVIDED THROUGH FILE *INPUT*

      CYCIT    1.0         MARCH     1995      EXECUTIVE ROUTINE FOR MARCHING IN TIME

      EOS      1.0         MARCH     1995      *LGMEOS3 * ... THERMOPHYSICAL PROPERTIES MODULE FOR WATER/THREE-COMPONENT AIR
      SAT      1.0      22 JANUARY   1990      STEAM TABLE EQUATION: SATURATION PRESSURE AS FUNCTION OF TEMPERATURE
      COWAT    1.0      22 JANUARY   1990      LIQUID WATER DENSITY AND INT. ENERGY AS FUNCTION OF TEMPERATURE AND PRESSURE
      SUPST    1.0      29 JANUARY   1990      VAPOR DENSITY AND INTERNAL ENERGY AS FUNCTION OF TEMPERATURE AND PRESSURE
      HENRY    1.0      25 FEBRUARY  1993      HENRYS LAW TO CALCULATE NCG BUBBLING PRESSURE OR DISSOLVED NCG MASS FRACTION
      VISW     1.0      22 JANUARY   1990      VISCOSITY OF LIQUID WATER AS FUNCTION OF TEMPERATURE AND PRESSURE
      VISCOSITY1.0      26 MARCH     1995      VISCOSITY OF GAS MIXTURE H2O, CO2, CH4 & N2, FOR GIVEN TEMPERATURE
      RELP     1.1         MARCH     1995      LIQUID AND GAS PHASE RELATIVE PERMEABILITIES AS FUNCTIONS OF SATURATION 
      PCAP     1.0       4 MARCH     1991      CAPILLARY PRESSURE AS FUNCTION OF SATURATION

      SHEAT    1.0      13   JUNE    1995      CALCULATES SPECIFIC HEAT Cv VALUES FOR GASES FOR GIVEN TEMP.

      BALLAN   1.0         AVRIL      1995      PERFORM SUMMARY BALANCES FOR VOLUME, MASS, AND ENERGY
      SECOND   1.1      17 AUGUST    1993      SUN SUBROUTINE FOR CPU TIMING FUNCTION
      TSTEP    1.0          MARCH    1995      ADJUST TIME STEPS TO COINCIDE WITH USER-DEFINED TARGET TIMES

      MULTI    1.2       4 OCTOBER   1993      ASSEMBLE ALL ACCUMULATION AND FLOW TERMS
      DIFF     1.0      22  MARCH    1995      DIFFUSION COEFFICIENTS FOR BINARY GAS MIXTURES AS FUNCTION OF TEMPERATURE
      LINEQ    1.0      22 JANUARY   1990      INTERFACE FOR THE LINEAR EQUATION SOLVER MA28
      MC19A                                    HARWELL SUBROUTINE FOR SCALING MATRIX

      LGMBAL   1.0      29 DECEMBER 1995      PRINT RESULTS FOR ELEMENTS, CONNECTIONS, AND SINKS/SOURCES
      CONVER   1.0       4 MARCH     1991      UPDATE PRIMARY VARIABLES AFTER CONVERGENCE IS ACHIEVED

      OUT      1.2      29 SEPTEMBER 1993      PRINT RESULTS FOR ELEMENTS, CONNECTIONS, AND SINKS/SOURCES
      WRIFI    1.0      22 JANUARY   1990      AT THE COMPLETION OF A TOUGH2 RUN, WRITE PRIMARY VARIABLES ON FILE *SAVE*
